February 1, 2008 4:19 WSPC/Trim Size: 9m x 6m for Proceedings 



conf2003 



ADVENTURES IN COULOMB GAUGE *t 



JEFF GREENSITE 

Physics and Astronomy Dept., San Francisco State University, 
San Francisco, CA 94117. E-mail: greensit@stars.sfsu.edu 



STEFAN OLEJNIK 

Institute of Physics, Slovak Academy of Sciences, 
SK-845 11 Bratislava, Slovakia. E-mail: fyziolej@savba.sk 



We study the phase structure of SU(2) gauge theories at zero and high temperature, 
with and without scalar matter fields, in terms of the symmetric/broken realization 
of the remnant gauge symmetry which exists after fixing to Coulomb gauge. The 
symmetric realization is associated with a linearly rising color Coulomb potential 
(which we compute numerically), and is a necessary but not sufficient condition 
for confinement. 



There are several reasons why Coulomb gauge may be interesting and/or 
useful in the study of the confining force. First of all there is the speculation 
by Gribov [1] and Zwanziger [2] that confinement in Coulomb gauge is due 
to instantaneous (dressed) one-giuon exchange. Secondly, the behavior of 
the color Coulomb potential, defined in Coulomb gauge, is an important 
element in the gluon-chain model of QCD string formation [3] . Finally, as 
we will see below, the confining property of the color Coulomb potential 
is associated with the unbroken realization of a remnant gauge symmetry, 
and this suggests a new order parameter for studying the phase structure 
of lattice gauge theories. 

We begin with the idea that confinement arises from one-gluon exchange 
in Coulomb gauge; specifically, from the instantaneous piece of the (AqAq) 
propagator 

(Ag(x)Ao(2/)} = P{x - y)J°''(5(xo - yo) + non-instantaneous (1) 
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where 



P{x-y)6 



ab 



1 

.V • D{A) 



1 

V-£»(A)_ 



(2) 



where Di{A) is the covariant derivative. This quantity is directly related 
to the Coulomb interaction energy in Coulomb gauge. 

Recall that the classical Hamiltonian in Coulomb gauge, H = Hgi^e + 
Hcouh has the form 

Hglue = ^ / (^*''" • + • -S") 

Hcoul = \l d^'xd^y p''{x)K'^\x,y;A)p\y) 



K'^\x,y;A) 



1 



1 



-I ah 



V-D{A)' ' 'V-D(A) 

P - Pmatter " 9 J ^k^k 



(3) 



Note that {K) is the instantaneous piece of the (^o^o) propagator. Gribov 
and Zwanziger argue that this propagator is enhanced by configurations 
at the Gribov horizon, defined as a boundary in function space where the 
operator V ■ D{A) aquires a zero eigenvalue. The conjecture is that this 
enhancement leads to a confining Coulomb potential, and therefore con- 
finement by one-gluon exchange. 

One objection to this idea is that it is difficult to see how the string-like 
properties of the QCD fiux tube, namely, the logarithmic growth of the flux 
tube cross-section (roughening), and the universal — 7r/12i? contribution to 
the static quark potential (the Liischer term), could arise from one-gluon 
exchange. On the other hand, as we will see below, the color Coulomb 
potential is an upper bound on the confining static quark potential. This 
means that confinement by one-gluon exchange is a necessary condition for 
confinement. 

Let 



(4) 



be a physical state in Coulomb gauge containing two static charges; \l/o is 
the ground state. Then 



AE = {^gg\H\^qg) - {^0\H\^0) 

= Vcoul (-R) + Ese 



(5) 



February 1, 2008 4:19 WSPC/Trim Size: 9in x 6in for Proceedings 



conf2003 



3 



is the expectation value of the excitation energy, where the i?-dependent 
Coulomb potential Vcoui (R) can only arise from the expectation value of the 
non-local pKp piece of the Hamiltonian. We want to address the following 
questions: Is Vcoui{R) confining? If so, is it asymptotically linear? If it is 
linear, then is the associated string tension acoui equal to string tension a 
of the static quark potential? Finally, we would like to study the effect, on 
the Coulomb string tension, of removing center vortices. 

We begin by defining the correlator, for SU(N) gauge theory in Coulomb 
gauge, of two timelikc Wilson lines 

G{R,T) = {j^TT[L\0,T)L{R,T)]) 

where 



L(f,T) = Pexp 



T 

dtA{x, t) 



(7) 



Note that L is a timelike Wilson line {not a Polyakov line) of time extent 
T. The existence of a transfer matrix implies 

G(i?,T) = ^|(*„|*,,)|^-^^"^ (8) 

n 

where the sum is over energy eigenstates, and A£'„ is the energy above the 
ground state. Denote 

V{R,T) = -^\og[G{R,T)] (9) 
Then its not hard to see that 

= Vcoul{R) + Ese 

= V{R,Q) (10) 
= \,mV{R,T) (11) 

1 — ^oo 

where AEmin is the minimum energy of the qq system, and V{R) is the 
static quark potential. The use of Wilson line correlators, in Coulomb 
gauge, to compute the static potential was first suggested by Marinari et 
al. [4]. 



while 
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With a lattice rcgularization, the self-energies Eg^ and E'^^ become neg- 
ligible, at large R, compared to the confining static potential V{R). Then, 
since AEmin < ^E, it follows that asymptotically 

V{R) < V,oui{R) (12) 

This inequality was originally derived by Zwanziger [5]. It implies that 
if confinement exists, it exists already at the level of dressed one-gluon 
exchange. 

With a lattice rcgularization, we have 



L(x,T) = Uoix,a)Uo{x,2a) ■ ■ ■ Uo{x,T) 
GiR,T) 



V{R,T) = ^log 



G{R,T + a) 



(13) 



so that 



lim V{R,0) = VcouiiR) + const. 

lim V{R,T) = V{R) + const. (14) 

where, in lattice units a = 1, 

V{R,0) = -log[G{R,l)] (15) 

Via lattice Monte Carlo we can then arrive at an estimate, exact in the 
continuum limit, of Vcoui{R) from V{R,0), and compare this to the static 
quark potential V{R). 

This procedure was carried out in ref. [6] for SU(2) lattice gauge theory. 
The result for V{R, 0) at /? = 2.5 is shown in Fig. 1 (upper set of data 
points). The two lines shown are best fits to the data by a linear, and by 
a linear + Liischer, potential. The data immediately answers three of the 
questions posed above: The Coulomb potential is confining, and it is linear 
(sec also [7]). However, it turns out that the associated string tension Ucoui 
is substantially greater than a, by almost a factor of three. This means that 
the QCD flux tube is not simply the static charges and their associated 
Coiilomb (longitudinal color electric) field. The minimum energy string 
state is more complicated than (4); it must also contain some constituent 
ghions, as in the gluon chain picture of string formation advocated in [3]. 

The same figure shows the effect, on V{R, 0), of center vortex removal. 
The center vortex theory of confinement has been studied very actively in 
recent years; the theory and the numerical evidence in its favor are reviewed 
in ref. [8]. Center vortices are identified by first fixing to an adjoint gauge, 
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DC 
> 



V(R,0)=-log[G(R.1)], [3=2.5 
with vortices 

0.483-11/(12 R) + 0.098 R 
0.392-1- 0.105 R 
without vortices 




Figure 1. Lattice approximation V{R, 0) (upper data set, open squares) to the color 
Coulomb potential Vcoui(^) + constant. Tiie lines are best fits to linear, and linear 
+ Liischer, potentials. The lower data set (open circles) for V{R, 0) is obtained after 
removing center vortices from the lattice configurations, as described in the text. 



and then projecting link variables to the center subgroup of the gauge 
group. An example is the direct maximal center gauge {— Landau gauge in 
the adjoint representation), where the procedure is to gauge fix to a local 
maximum of 



(16) 



and then to project each link to the closest center element, e.g. for SU(2) 



Ufj,(x) Zf_,{x) = signTr[C/p(x)] 



(17) 



Vortices arc removed from a given lattice configuration by multiplying the 
adjoint gauge-fixed configuration by the projected configuration, i.e. 



(18) 



In ref. [9] it was shown that after vortex removal the string tension van- 
ishes, chiral symmetry breaking is eliminated, and each vortex-removed 
configuration has zero topological charge. 

One can then ask what vortex removal does to the color Coulomb po- 
tential. In this connection, a relevant fact is that thin center vortices can 
be shown [10] to lie on the Gribov horizon, which is thought to play an im- 
portant role in the enhancement of the Coulomb energy. In our numerical 
study, the modified configuration U'^{x) is gauge- fixed to Coulomb gauge, 
the timelike link correlators are calculated, and V{R, 0) is extracted. The 
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result a,t (3 — 2.5 is shown in the lower data set (open circles) in Fig. 1. It 
is clear that vortex removal completely removes the confining property of 
the Coulomb potential. Further details can be found in ref. [6]. 

Of course, pure gauge theory at zero temperature is a special, albeit 
very important, case. In general gauge theories can exist in various phases, 
and it is useful to characterize these in terms of the distribution of the 
electric field emanating from a static isolated source. 

• Massless Phase: The electric field is spherically symmetric, and 
falls off like 1/R'^. This is the case for compact QED4, and for 
lattice SU(N) gauge theory in £> > 4 dimensions, at weak couplings. 

• ConGncd Phase: The color electric field is collimated into a 
flux tube; global center symmetry is unbroken. Examples include 
SU(N) pure gauge theories at low temperature, and SU(N) gauge 
theories with matter in the adjoint representation of the gauge 
group. 

• Screened Phases: There is a Yukawa-like falloff of the color electric 
field. This happens in SU(N) gauge theories when the Zn center 
symmetry is broken spontaneously, as in high temperature gauge 
theory and gauge-Higgs theories with the Higgs field in the adjoint 
representation. Gauge theories with only a trivial center symmetry 
(consisting only of the identity element) are also in the screened 
phase; these theories include SU(N) gauge theories with matter 
in the fundamental representation, and G2 gauge theory with or 
without matter fields. 

For the purpose of studying Coulomb energy in these various phases, we 
find it useful to introduce a new order parameter, related to the realization 
of a remnant gauge symmetry in Coulomb gauge. 

Let C/p(x) be a lattice configuration fixed to Coulomb gauge. Note that 
the Coulomb gauge condition is preserved by the gauge transformation 

Uk{x,t) ^ g{t)Uk{x,t)g^{t) 

Uo{x,t) ^ g{t)Uo{x,t)g^{t+l) (19) 

On any time slice, this is a global transformation, and therefore can be 
spontaneously broken in the following sense: At any fixed time t, in the 
infinite volume limit, the average of timelike link variables Uo{x,t) is non- 
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zero in any thermalized configuration. This means that 
lim G{R, 1) > 
lim V{R, 0) = finite const. 

-R— >oo 

cTcoul = (20) 

in the broken phase. Therefore Coulomb confinement or non-confinement 
can be understood as the symmetric or spontaneously broken realization, 
respectively, of the remnant gauge symmetry in Coulomb gauge. 

We now introduce the order parameter Q, as the modulus of the spatial 
average of timelike links, i.e. 

Q=(^^Tr[US^{t)U^^\t)]j (21) 

On general grounds 

Q — c+ ^ 'th I ~ ^ symmetric phase 
\ c > broken phase 

Thus Q > in the infinite volume limit implies that Vcoui{R) is non- 
confining. Since Vcoui{R) is an upper bound on V{R), this implies that 
Q = is a necessary (but not sufficient) condition for confinement. 

It is useful to try out this order parameter in compact QED4, where 
there is a transition from the confining to the masslcss phase at /3 « 1. Fig- 
ure 2 shows our results for Q vs. the root inverse 3- volume L~^^^ at /? = 0.7 
(confining phase) and /? = 1.3 (massless phase). In this case the Q param- 
eter seems to nicely distinguish between the two phases, extrapolating to 
zero only in the confined phase. 

Next, we consider SU(2) gauge-Higgs theory with a "frozen" Higgs field 
in the adjoint representation. The lattice Lagrangian is 

s = pY^^Ti[uuu''u''] 

plaq 

+ i E -^^ W-/-'!^ + li)Tr[a'^U,{x)a''Ul{x)] (23) 

where is a real 3-component field satisfying the constraint X^^aC*^")^ ~ 
This is a theory with a confining, center symmetric phase, and a non- 
confining phase with spontaneously broken center symmetry. Our finding 
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Figure 2. Extrapolation of Q to infinite volume in QED^, for /9 = 0.7 (confining phase) 
and /3 = 1.3 (massless phase). 



is that the transition hne in the /? — 7 phase diagram corresponding to the 
remnant symmetry-breaking transition is identical to the transition line for 
confinement-deconfinement, mapped out long ago by Brower et al. [11] from 
measurements of the plaquette energy. In the confined phase we find Q = 
(when extrapolated to infinite volume), and Q > in the Higgs phase, as 
indicated schematically in Figure 3. 





Figure 3. Phase diagram of the SU(2) adjoint Higgs model. The plaquette energy and 
the remnant symmetry order parameter Q locate the same transition line between the 
confined and Higgs phases. 



One might guess that the transition from the confined to the deconfined 
phase is always accompanied by remnant symmetry breaking. Surprisingly, 
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this turns out not to be true. We have also computed V{R, 0) and Q in 

the deconfined phase of pure SU(2) lattice gauge theory, with the results 
shown in Figs. 4 and 5. This data was taken at (3 = 2.3 on lattices with 
time extension of two lattice spacings, well within the deconfined phase. 
Yet the Coulomb potential is clearly linear and confining at large lattice 
volume, while the extrapolation of Q to infinite volume seems compatible 
with zero. A possible reason for this behavior is the fact that K{x,y;A), 
whose expectation value gives the instantaneous Coulomb propagator, de- 
pends only on the spacclike components Ak at a fixed time. On the lattice, 
this translates to dependence only on spacelike links on a time slice. But 
we know that spacelike links on a time slice are a confining ensemble even 
in the deconfined phase, since spacelike Wilson loops are known to have an 
area law falloff at any temperature. If the Coulomb propagator depends 
only on the confining properties of spacelike links, then it is not so sur- 
prising that the Coulomb potential is confining in the deconfined regime 
(nor is this a paradox: the Coulomb potential is only an upper limit on 
the static potential). A test of this explanation is to remove the confining 
properties of the spacelike links by removing center vortices, via the de For- 
crand/D'Elia procedure explained above. Then one expects the Coulomb 
potential to be non-confining, and this is. in fact, what is observed. 

Finally, we study a gauge-Higgs system with the radially frozen Higgs 
field in the fundamental representation. For the SU(2) gauge group, the 
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Figure 5. The Q parameter vs. root inverse 3-volume in the deconfined phase, P = 2.3 
and Lt = 2 lattice spacings. 



lattice Lagrangian can be expressed as [12] 
S = (3^^Ti[UUU'^U^] 

plaq 

+ 7^iTV[0t(a;)C/^(x)<^(x + M)] (24) 

with (j) an SU(2) group-valued field. This is a theory with only a screened 
phase; it can be proven that no transition to a confined phase is possible 
[13]. There is a first-order phase transition line in the /3 — 7 phase dia- 
gram, but this line has an endpoint, and does not divide the diagram into 
thermodynamically separate phases. 

The remnant symmetry transition line coincides with the (thermody- 
namic) line of first-order transitions found by Lang et al. [12], but it then 
extends beyond the thermodynamic line all the way to /3 = and 7 = 2. 
This line divides the phase diagram into Q = and Q > regions, as indi- 
cated schematically in Fig. 6. In Fig. 7 we plot Q vs. 7 at /? = 0. If Q were 
the magnetization of an Ising spin system, this would surely bo a second 
order phase transition, with the solid line in the figure representing the 
infinite volume limit. Nevertheless, there is no thermodynamic transition. 
At /3 = one can easily compute the bee energy exactly, which is found to 
be 

(25) 



F(7) =41/log 



2/1(7) 
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Figure 6. Phase diagram of the fundamental Higgs model. There is a thermodynamic 
transition and a Q transition along the solid line, but a non-thermodynamic transition 
(Kertesz line) in Q along the dashed line. 



SU(2) with fund. Higgs, M 
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0.6 

0.4 
0.2 



1 2 3 4 5 6 7 



Figure 7. Q vs. 7 at /3 = in the SU(2) fundamental Higgs model, on 8* and 16* 
lattices. The solid line is the presumed extrapolation to infinite volume. 



This expression is perfectly analytic at all 7 > 0. On the other hand, 
a strong-coupling analysis of G{R, 1) at fixed (3^1 [10] arrives at an 
exponential decay to zero as i? ^ 00 at small 7, but a non-zero large- 
R limit at large 7. This implies a symmetry-breaking transition at some 
critical value 7cr(/3), which motivated our numerical study of Q in this 
model. 

The remnant symmetry breaking transition in the gaugc-Higgs system, 
in the absence of a thermodynamic transition, is probably an example of a 
Kertesz line [14] in statistical mechanics. This possibility was first suggested 
by Langfeld [15], who discovered remnant gauge symmetry breaking in 
Landau gauge in a closely related model (see also Satz [16]). Kertesz lines 
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are associated with percolation transitions, and it is natural to ask, in 
this case, what is percolating. Recent investigations [17] indicate that in 
the gauge-Higgs system, the Kertesz line locates center vortex percolation 
transitions. 

This concludes a summary of our investigations of phase structure in 
lattice gauge theory, as seen by Coulomb energy and remnant symmetry. 
Our study has uncovered Coulomb "over-confinement" {acoui ~ 3(t) in the 
low temperature confined phase, the persistence of Coulomb confinement 
in the deconfined phase, connections between vortex and Coulomb confine- 
ment, and (in accord with Langfeld [15]) symmetry breaking in the absence 
of a thermodynamic phase transition. These aspects of non-perturbative 
gauge theory are somewhat surprising, and merit further study. 
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